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I ■ Abstract. We present a class of cooperative sequential adsorption models on a 

^ , Cayley tree with constant and variable attachment rates and their possible applications 

for ionic self-assembly of thin films and drug encapsulation of nanoparticles. Using 
the empty interval method, and generalizing results known from reaction-diffusion 
^ [ processes on Cayley trees, we calculate a variety of quantities such as time-dependent 

surface coverage and time-dependent probabilities of certain particle configurations. 
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1. Introduction 

Nanotechnology research is a highly interdisciphnary area of study that rehes on 
fundamental understandings of physics, chemistry, biology and materials science. Layer- 
by-Layer self-assembly of nanoparticles is an innovative, widely used technique in 
nanotechnology studies with applications in microelectronics, optical coatings and 
biology [1]. In this context, non-equilibrium statistical physics methods and models 
can shed some light on the cooperative behavior of multi-particle systems. 

One-dimensional sequential adsorption models have been studied thoroughly in 
different physics contexts [2], [3]. The two basic one-dimensional models are known 
in literature as random sequential adsorption (RSA), when the adsorption sites are 
chosen randomly, and cooperative sequential adsorption (CSA), when the adsorption 
mechanism is influenced by the local environment. Despite numerous studies of one- 
dimensional models, adsorption in two dimensions is not as well understood. There are 
quite a few computational papers [2] on the matter, but very few analytical solutions 
exist for the general two-dimensional case. 

The adsorption of particles is exactly solvable in higher dimensions only for a finite 
tree-like lattice called a Cayley tree. A Cayley tree is a connected, cycle-free graph with 
each node connected to z neighbors, where z is called the coordination number (Fig. 1). 
For z = 2, the tree reduces to the one-dimensional line. It has been proven 

that a Cayley tree structure with z = 4 well approximates a regular two-dimensional 
lattice for certain cases, such as monomer and dimer deposition. For example, the Cayley 
tree predicts an 88.9% dimer coverage while the regular lattice predicts a 90.8% dimer 
coverage (result obtained from computer simulations [2]). Recently, analytical results 
were reported for the random sequential process [5] and reaction-diffusion processes on 
Cayley trees and Bethe lattices [6], [7], [8]. 

In this paper, we study a cooperative sequential model of monomer and dimer 
deposition with both constant and variable attachment rates on a Cayley tree. Two 
experimental topics motivate our paper. One is the understanding of the self-assembly 
mechanism of charged nanoparticles on a glass substrate [9j. Known in literature as 
ionic self-assembled monolayers (ISAM), this technique has been used successfully in 
the making of antireflective coatings [lOJ. The properties of these coatings depend 
on the surface coverage of the substrate. During the manufacturing process, it is 
highly desirable to know the analytical relationship between the index of refraction 
and the particle density of the surface. We study the ionic self-assembly deposition 
from the point of view of stochastic cooperative sequential adsorption models on a 
Cayley tree lattice. For the purpose of building a theoretical model for this process, 
we emphasize some of its essential elements: i) The deposition process is stochastic, so 
a random sequential adsorption model is justified; ii) The nanoparticles deposited are 
electrically charged, so a cooperative sequential adsorption model with rates dependent 
on the nearest neighbor occupation is appropriate; iii) The surface on which particles 
are being deposited is covered with a polyelectrolyte. We consider a Cayley tree 
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as an approximation for this surface. This choice is particularly appropriate when 
the polymer used is poly(amidoamine) (PAMAM), which has by design a tree-like 
structure. For different polyelectrolytes, such as poly(allylamine hydrochloride) (PAH) 
or poly(diallyldimethylammonium chloride) (PDDA), this approximation works well 
enough for special cases that will be addressed in section 3. 

The other experimental interest involves the versatile properties of synthetic 
polymers called dendrimers, with potential use as a novel drug delivery mechanism via 
drug encapsulation [II],[I2]- Dendrimers are perfect physical examples of Cayley tree 
structures. They are highly branched, spherical polymers that consist of hydrocarbon 
chains with variable functional groups attached to a central core molecule. Due to 
the precise control that can be exerted over their size, molecular architecture, and 
chemical properties, dendrimers have great potential in the pharmaceutical industry as 
effective carriers for drug molecules. These new synthetic polymers are able to carry both 
targeting molecules and drug molecules to cancerous tumors, minimizing the negative 
side effects of medications on healthy cells. In order to describe mathematically the 
encapsulation process of drug molecules, we propose a stochastic deposition model with 
variable attachment rates. 

Our paper is structured as follows. In section 2 we define our general model for 
both constant and variable rates. In section 3 we discuss in more depth the ionic self- 
assembly mechanism and calculate the quantities of interest, such as the coverage of the 
surface in the steady state, the probability of having clusters of connected empty sites 
of a certain size, and the probability of having a certain shell particle distribution. In 
section 4 we address the drug encapsulation model and present a set of analytical results 
for the case of variable attachment rates. We conclude in section 5 with a summary of 
our results and some open questions. 

2. Model definition 

Random sequential adsorption of a mixture of monomers and dimers on a Bethe lattice 
was solved by Cadilhe and Privman [5]. Using the empty interval method, they 
found exact analytical expressions for the time evolution of clusters of connected empty 
sites. We address the question of adsorption of monomers and dimers with cooperative 
effects due to different neighbor occupation. We consider the following general processes 
defined on a Cayley tree: 

• o — > • • with rate ri 
o o — > • o with rate r2 
o o — )■ • • with rate rs 

where "•" are filled sites and "o" are empty sites. The same processes, with constant, 
positive rates, were considered in the context of diffusion-reaction problems by Matin 
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et a/.[n] and Alimohammadi et al. [7j. 

In the following sections, we consider two distinct cases for the deposition rates ri, 
r2 and rs: i) constant rates, to model the ionic self-assembly process; ii) variable rates, 
to model the drug encapsulation mechanism using dendrimers. 

For both cases, some general information regarding Cayley trees is needed. 
According to the mathematical definition, a Cayley tree of order z is defined in the 
following way [13]. Given a root vertex 0, we link this vertex to z new vertices by means 
of z edges. Shell £ = 1 is made of the first set of z new vertices. Each vertex on shell 
£ = 1 is linked to z — 1 new vertices, and so on. Fig. 1 shows the first 3 shells for a 
Cayley tree of coordination number z = 4. Cayley trees are finite trees with boundaries 
defined by the last shell; in the case presented in Fig. 1 this is shell L = 3. The number 
of vertices on the shell is Hi = z{z — lY~^, and the total number of vertices of the 
Cayley tree with L shells is = z{{z - 1)^ - l)/{z - 2). 

We calculate the following quantities: 

• En(t), the probability of finding a cluster of n-connected empty sites inside the tree, 
regardless of the occupation of the rest of the nodes. 

• PL{m;t), the probability of having m particles on the last shell L, with all other 
shells empty. 

• Pi{m;t), the probability of having m particles on the shell i < L, with all other 
shells empty. 

• P{mi] 7712; ■■'mi] mi] t), the probability of having {mi; m2; ■■me} particles distributed 
on interior shells and particles on the last shell, provided that there is one shell 
empty between filled levels. 

3. Simple model of ionic self-assembly: cooperative sequential adsorption 
with constant rates 

Ionic self-assembly is a relatively new technique ^14j, [I5j that allows detailed structural 
control of materials at the molecular level, combined with ease of manufacturing and low 
cost. The ISAM process allows one to deposit alternating layers of cations and anions by 
dipping the substrate in aqueous solutions of the appropriate ions, as illustrated in Fig. 
2. Because it is a dipping process, any exposed surface is homogeneously coated, allowing 
highly uniform, conformal coatings on irregular shapes. The basic building block for 
the film is a cation/anion bilayer, which may consist of either two polyelectrolytes (a 
polycation and a polyanion), or a polyelectrolyte and a nanoparticle, or two different 
nanoparticles. The thickness of a bilayer is a function of the diameter of the nanoparticle 
and the packing of the particles from layer to layer. Fig. 3 shows a Si02/PDDA 
bilayer made with spherical silica particles with nominal sizes of approximately 45 nm, 
self-assembled on a glass slide. Some dimers are visible as well. 

We model the ISAM mechanism as a cooperative sequential adsorption of 
nanoparticles with constant deposition rates ri, r2 and r^ on a Cayley tree. We 
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performed Monte Carlo simulations for both a Cayley tree of coordination number 
z = 4 and a fully two-dimensional lattice and monitored the filling process of both 
lattices. The simulation results show that for similar values of ri and r2, there is a 
very good match between a Cayley tree and a fully two-dimensional lattice. The results 
start to differ once ri and r2 are significantly different. The attachment rates can be 
picked to fit the experimental data. A pH = 9, for example, translates into a higher 
attachment rate for monomers than pH = 8. This model was discussed in the context 
of reaction-diffusion reactions on a Cayley tree by Matin et al. ^ and Alimohammadi 
et al. [7]. 

3.1. Calculation of probability distribution of clusters of n- connected empty sites using 
the empty interval method 

For the ionic self-assembly process, one of the main concerns is the time dependent 
surface coverage, and its properties in the steady state. Let -P(On) = En be the 
probability of finding a cluster of n connected empty nodes, regardless of the state of 
the bordering sites. For clarity, in Fig. 4(a) we present a cluster of four connected 
empty sites, E^. From now on, we assume that an "n-cluster" refers to a cluster of n 
connected sites, unless specified otherwise. Let P(» — On) be the probability to have 
an n-cluster that is empty connected to an occupied node via one single link. In this 
case, we have the following conservation of probability equation [16j : 

- On) + P(o - On) = P{On) (l) 

which leads to 

Pi* - On) =En- En+1 (2) 
The equation that governs the probability of having an ra-cluster of empty sites is: 

^ = -{niz - 2) + 2)(riP(. - On) + (^2 + r3)K+i) - (2r2 + r,){n - 1)E„ 
dt 

(3) 

which can be simplified to: 

dE. 



dt 



- = -[n{n{z - 2) + 2) + (2r2 + r3)(n - l)]En - (r2 + r^ - ri){n{z - 2) + 2)E„ 



For the special case of ri = r2 + r^, and a completely empty lattice at the initial time, 
the solution for En is: 

For the general case of ri 7^ r2 + r^, using the standard method presented in [T7] . 
we obtain the same results as Matin et al. |6j. We use the following exponential ansatz: 

En{t) = E,{tmt)r-' (6) 
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We substitute this into Eq. 4, and we get a system of coupled differential equations 
by separating the initial equation in terms linear in n and independent of n: 

'^'^^^^ = - m) - f^a^itf (7) 



dt 

dEi{t) _ 
dt ~ 

with the notations: 



[zr, + (8) 



I3 = {z- 2)ri + 2r2 + rg (9) 

(z-2)(r2 + r3 -ri) 

a = ^ (10) 

{z - 2)ri + 2r2 + rg 

With the appropriate initial conditions {En{t = 0) = 1, empty lattice), we recover 
the following results from [6]: 



1 + a{l — e P*) 

Enit) = e-^^-("-i)/^^( 1 -,n )^^"-' (12) 

1 + a(l — e P^) 

Fig. 4(b) displays the probability of finding a cluster of four empty, connected 
nodes (n = 4) for a Cayley tree with z = 3, and arbitrary deposition rates ri = 0.2, 
r2 = 0.4, r3 = 0.1. The initial lattice is empty, and the probability approaches zero in 
the steady state. 

We can also answer the question of the particle density of the jammed state. When 
all three attachment rates are non-zero, the final state will be completely full. For the 
case of ri = 0, for example, the final state is 

^ / (r2 + rs) {z - 2) e-(^-^+-3)i + ^3 _ ^ (^^ + ^3) \ 

^ ' V -2r2-r3 J ^ ' 

In Fig. 5(a) we present the time dependence of the particle density for 2 = 4, 
Ti = 0, r2 = 1, and r^ = 0.1. The final coverage is 76%. In Fig. 5(b) we have the 
steady state coverage {t — )■ 00) as a function of coordination number z for the same set 
of parameters. The coverage increases with time, and tends asymptotically to 100% as 
z grows larger. 

3.2. Probability distributions of shell occupation 

For the ionic self-assembly mechanism, it is interesting to know not only the final surface 
coverage, but also the probability of having certain particle patterns. The results 
presented below are apphcable for the case of a "snowflake" type polymer, such as 
PAMAM. We assume below an ideal case of particle deposition on only one "snowflake" 
of generation L. Alimohammadi et al. [7] addressed the possible reaction-diffusion 
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processes on a Cayley tree, and showed that the only reactions amenable to exact 
solutions are the ones listed in our model. Below, we discuss some of their results in the 
context of particle deposition with constant rates, and we also present a more general 
solution to the problem. For convenience, we remind the reader of the quantities of 
interest: 

• Pi{m;t), the probability of having m particles on the last shell L, with all other 
shells empty. 

• Pi{m\t), the probability of having m particles on the shell i < L, with all other 
shells empty. 

• P{mi] m2] ..mi] rriL] t), the probability of having {mi; 7712] ..m^} particles distributed 
on interior shells and particles on the last shell, provided that there is one shell 
empty between filled levels. 

The empty shell between filled consecutive shells is necessary in order to ensure a closed 
set of differential equations [7j . 



a) Internal shells i < L 



The master equation that governs particle deposition is: 
d PAui' i) 

— = {ni—m+l)zr2Pi{m—l]t) — \mzri + {2r2+r3){N—rnz)]P£{rn]t){lA) 

dt 

with n£ being the number of sites on the i^^ shell and N = Y,i'=i f^e the total 
number of sites. 

[ ziz-lY-^ for t>l 
"^=|l for £ = 

The first term on the right hand side of Eq. 14 accounts for the processes that lead 
to configuration Pe{m; t). Starting with m — 1 particles on shell i, there are {ne — m + l)z 
possibilities for a particle to be deposited on shell i. The negative terms in the master 
equation account for the disappearance of configuration Piijn] t) via particle depositions 
on shell I or on the other interior shells. 

For simplicity, we introduce the following notations: 

a = ri — 2r2 — (16) 
b = N{2r2 + r3) (17) 

With these notations, Eq. 14 can be rearranged as: 
P (tti' t) 

— ^^-^ = {ni-m+ l)zr2Pe{m - 1; t) - (mza + b)Pi{m; t) (18) 

The method used is recursive, similar to [7]. The starting point is the case of m = 0, 
when the equation becomes: 

dPAO;t) , , , 

^^ = -6P,(0;t) (19) 
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We assume an empty initial configuration, so P^(0; 0) = 1. The solution is: 

P,(0;t)=e-^* (20) 

The solution for Pf(0;t) becomes the "starter" for the equation for P^(l;t) and so on, 
to get the following general solutions: 

^ ^."(-irr(^-1.)e--(l-e-r „ ^ „ 

i (— n^) ml a"* 

P,(,„;0='''-''-'"'-V:rt'"-'")^-"fo.a = (22) 

m\V[-n() 

Fig. 6(a) shows the time dependence of the probability of having four particles in 
the third shell of five, with all other shells empty, for z = 4, vi = 0, r2 = 1, and r^ = 0.1. 
The probability reaches a maximum of roughly 7.08 x 10~^ at i = 4 x 10~^ in arbitrary 
units. 

Larger r-values shrink the time scale and condense the curve. Smaller r-values 
elongate the time scale, that is, the lower the rate of deposition, the longer it will take 
to get to the desired number of particles in the shell. Of the three, r2 and have the 
greatest influence over the curve, altering both the time scale and the order of magnitude 
of the probability for the particle distribution in the shell. We can also find analytically 
tmax, the peak of the distribution: 

((^-1)^-1)^ 



In 



za 

b) Last shell L 

The master equation that governs these processes is (set ^ = 1 in Eq.l4): 

dPL{m;t) 



dt 

with the solutions: 



{nL-m + l)r2PL{m -l;t)- {ma + b)PL{m; t) (24) 



r (— n^) ml a^ 

, (r2t)"(-l)™r(m- nL)e-^* , , 

Pl m; t) = > ^ / for a = 26 

m!r (-n^) 

Fig. 6(b) shows the time dependence of the probability of having four particles in 
the final (L = 5) shell with all other shells empty for z = 4, ri = 0, r2 = 1, and — 0.1. 
The probability reaches a maximum of roughly 2.01 x 10~^ aX t — 3.95 x 10~^. Notice 
the overall similarity of this curve to the curve for interior shells. The relationships 
between the r-values and the curve are the same as in Fig. 6(a). Note that the order 
of magnitude for the probability here is 10~^ or less, compared to 10^"'' for an interior 
shell. This is the result of the outer shell always having the most possible spaces to 
occupy, so it is more likely that it will fill to a certain point with the rest of the tree 
empty than for any interior shell to fill to that same point under the same conditions. 
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As before, we can find an analytical expression for the time when the maximum of the 
distribution is reached: 

a 

3.3. General solution for P{m,t) 

Consider a case of m = mi + m2 + rn^ + .. particles being distributed on interior shells 
with at least one shell empty between them. We want to calculate the probability of 
this particle distribution function of time. 

The general form for the equation for one shell occupied with all others empty can 
be expressed as: 

^^^^^ = r2z{ni -mi + l)P(m£ - 1) - [mza + h)P{mi), (28) 
for any inner shell £. 

Using Eq. [28] as an example, we move the h term from the right to the left side of the 
equal sign, allowing us to rewrite this equation as: 

+ = exp(-6t)^^^^^PM^^ (29) 

= r2z{ni — me -\- l)P{me — 1) — mzaP{mg). 

Noting that the time dependence of the totally unoccupied structure is: 

P(0) = Pi(0) = P(0, 0, ■ ■ ■) = exp(-6t), 

and defining the time-dependent functions Q{me) = P{m£)/ P{0) , we can rewrite 
Eqns. 25 and [281 in terms of the Q's as: 

^^^^^ ^ 'Z'?^^ ^ 

— = r2z{n£ — m£ + l)Q{me — 1) — mzaQ{me) (30) 

at 

For the case of an outer shell, the equation is: 
dQ(mij^ 

^ = r2{nL -mL + l)Q{mL - 1) - maQlmi) (31) 

Now consider the case of three occupied shells - two inner shells and the outer shell 
separated by at least one totally empty shell - and examine the equation satisfied by 
the product Q{mi,mk,mL) = Q{mi)Q{mk)QL{mL)- 

ot ot ot 
— Q[mi)QL[mk) 

= r2z{ni + 1 - mi)Q{mi - 1, m^, m^) + r2z{nk + 1 - mk)Q{mi,mk - 1, mj^) 
+r2{nL + 1 - mL)Q{mi, m^, - 1) - {az{mi + m^) + zmL)Q{me, mk, mi) 

(32) 
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This is completely equivalent to the master equation written directly for three occu- 
pied shells upon making the substitution Q{mi,mk,mL) = P{mi)P{mk)Pi{mL)/P{0)^. 
This procedure is trivially extended to cover any number of occupied shells, thus gen- 
eralizing the following result: 

o/ +^ Ili€{l,2,..£}P{'rni;t) 

P{mi;m2;..me;mL;t) = — — ■ — (33) 

P(0; 0; t)^'"*"' 

In the set {1, 2, ..£}, "1" is the first occupied shell, "2" is the second occupied shell, 
etc., with the assumption that there is always at least one empty shell between occupied 
shells. Lfotai is the total number of occupied levels (including the last external shell). 

4. Nanoparticle deposition with variable attachment rates 

The shell- dependence of the deposition rates is particularly relevant for the drug 
encapsulation mechanism via dendrimers. There are two methods of drug delivery 
using dendrimers: attachment of drug molecules to the outer functional groups and 
the encapsulation of drugs within the cavities created by separate branches. For 
example, dendrimer surface chains can form covalent bonded complexes with anti- 
cancer molecules such as cisplatin [TT]. The release of these molecules can be controlled 
by manipulating the rate at which these dendrimer-drug bonds are degraded. The 
potential load of each dendrimer carrier can be easily varied by adjusting the branch 
multiplicity of the dendrimer generation. Another mechanism of dendrimers to load 
molecules is through molecular encapsulation. The branches of the dendrimers form a 
dendritic box around the encapsulated molecule, which can protect sensitive molecules 
from unfavorable physiological environments, and vice-versa. 

The encapsulation process is complex, and depends on the type of dendrimer used, 
the type of drug molecules encapsulated, and other physical or chemical factors. We 
propose a minimalist model that considers the drug molecules as generic monomers 
capable of attaching to the available nodes. To capture the fact that it is less likely to 
have drug molecules present in the outer shells, we define shell-dependent attachment 
rates. We treat the surface of the dendrimer separately, since all nodes are available 
for attachment. Fig. 7 presents a sketch of the encapsulation of drug molecules and 
attachement of drug molecules onto the surface of a dendrimer. Dendrimers can carry 
different types of nanoparticles (illustrated by the different shapes in particles in the 
picture). 

In order to model the encapsulation mechanism of drug molecules much larger than 
the dendrimer nodes, we consider here a special case where the deposition of monomers 
is forbidden if the neighboring connected site is occupied. Dimer deposition is not 
being considered. The long time behavior of such a model is more interesting because 
eventually the system will reach a "jammed" state with a certain percentage of empty 
sites. Because of their size and their electrostatic interactions, drug molecules cannot 
be in close proximity, and this fact is reflected in our model by obstructing occupation 
of adjacent connected sites. Also, inside the dendrimer tree, the closer to the core the 
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drug molecules are, the more available space they have, and the more likely they are 
to be there. For our model, this translates into shell-dependent deposition rates that 
decrease for shells further from the core. For the drug encapsulation model, what is 
physically relevant is the particle distribution on the shells. We use the same analytical 
method as the one for the constant rates case. 

4.1. Probability distributions of shell occupation 

a) Monomer deposition on the last shell L: variable rates 

The following master equation describes the monomer deposition on the last shell 

L. 

^El^^ = r(L)K -m + !)PL{m - - [r(L)K - m) (34) 

L-1 

+ E^^'(^(0 + K^'-l))]^L(m;t) 
i'=i 

with ni being the number of sites on the shell. 

r z{z - for l>l 

= for / = ^^'^ 

The first term on the right hand side accounts for the processes that lead to configuration 
Pf,(m; t). Starting with m — 1 particles on layer L, there are (n^ — m-|- 1) possibilities for 
a particle to be deposited on shell L. The negative term in the master equation accounts 
for the disappearance of configuration PL{m; t) via particle depositions on shell L or on 
the interior shells I — 1...L — 1. 

We assume that r{L) = 1, and for now we consider for the rest of the shells 
r(/) = A', with A < 1 as a function of the shell number. The choice of r(/) is consistent 
with electrostatic screening and with the distribution of drug molecules encapsulated in 
the dendrimer. 

L-1 

l'=l 

The sum on the right hand side can be written as: 

j:n,ir{l') + r{l' - 1)) = z{\ + 1) Y.iK^ - ^)Y~' (37) 
l'=l v=i 

If we work under the assumption that A < then the sum in Eq. 37 is the 
partial sum of a geometric series: 
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<A + 1) UXiz - 1))''- = .(A + 1) V^^T^'" 



) 

(38) 



Eq. 36 can be written in a more compact form, 



dPL(m:t) , , ... , (A + l)(-2 - n^A^"^), „ , , 
^ ^ = (nL-m+l)PL(m-l;i)-[(ni-m)+^ ^ -\PL{m-t){m) 

The general solution of this equation is of the form: 

oo 

PL{m-t)= J2b^Pre'-' (40) 

m=0 

The eigenvalues, given by the coefficients of PL{m;t) are: 

= -[{riL - m) + ^ _\^^ _ -^^ ] (41) 

For the assumptions made, these eigenvalues are always negative, so the probability 
of having m particles on shell L while the other shells are empty decreases with time, 
as one might expect. Given initial conditions, and set system size, we can find the 
corresponding eigenvectors P£". To achieve the same goal, we use a recursive method. 
For convenience, we use the following notations from now on: 

= z{z - 1)^-^ (42) 
, . {X + l)(z-ziz-l)''-'X'^-') 

7 = . (- ir- + \-A(.-i) — - 



Eq. 39 then becomes: 
dPL{m;t) 



= (riL -m + l)PL(m - l;t) - (7 -m)PL(m;t) (44) 
at 

First, we start with m — 0: 

|Pi(0;t) = -7PL(0;t) (45) 
We assume an empty initial configuration, so Pl{0;0) — 1. The solution is: 

PL(0;i)=e-^* (46) 

With this solution, we go into the equation for PL{l;t) and solve it to get 

Pi(l;i)=nz,(l-e-*)e-(^-^)* (47) 
Following this method, we find a general form for the probability as: 

P(n^-f)- r(m-n,) (e-*-l)-e-(— 

r(-n.)m! ^^^^ 
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Fig. 8(b) displays the probability distribution of the deposition of two drug parti- 
cles (m = 2) on the final shell of a five generation (L = 5) Cayley tree with 2; = 3 and 
arbitrary inner-shell-dependent deposition rate of r{l) = 0.4'. Note that r{L) = 1. The 
probabilities here are two orders of magnitude higher than those of the third distribution 
for the same m (see Fig. 8(a)). The probability reaches a maximum of roughly 0.17 at 
t — 0.037, again in arbitrary units. 

b) Monomer deposition on interior shells I < L: vairiable rates 

The same steps can be followed to solve for Pi{m;t), the probability of having m 
particles on shell / < L, provided that the rest of the shells are empty. In this case, the 
master equation that needs to be solved is: 

^^^^ = ril)z{ni -m + l)Pi{m -l;t)- [ul - mzr{l) (49) 

+ Y.Mr{l')+r{l' -l))]Pi{m-t) 

i'=i 

with ni being the number of sites on the shell. 

The equation can be written in a more compact form: 



dPi{m;t) 



dt 



= Xz{ni -m + l)Pi{m -l;t)- [ul - mzX (50) 
+ l-A(.-l) 

The first term on the right hand side represents the process of adding a particle in the 
I shell. There are (n^ -|- 1 — m) empty sites in the I shell before the addition. Each 
such site is part of z pairs, and the rate at which particles attached to /-shell sites is 
r(/). Thus, the term has the form +{ni + 1 — m)zr{l)Pi{m — 1). Now, adding an L-shell 
particle when there are no particles there can happen in any of the Ul empty sites, 
each of which is part of only one pair. The filling rate is 1, producing the next term 
—riLPiim). Each inner shell /' (except shell /) is totally empty, thus has riii tillable sites, 
each a part of z pairs, filled at rate r(/'), producing —niizr{l')Pi{m). For shell I this rate 
is —{nil — m)zr{l')Pi{m). 

In terms of 7, the equation can be rewritten as: 

^El^ ^ (^^ + 1 _ m)zX^Pi{m - 1) - (7 - mzX^)Pi{m). (51) 

The matrix associated with the master equation is tridiagonal, and the eigenvalues 
are the coefficients of Pi{m; t): 

em = -{l-mz\') (52) 
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Using the same method as in the previous section, we find for P;(m; t) the following 
solution: 

with ni = z{z — iy~^. 

Fig. 8(a) displays the probability of deposition of two drug particles (m = 2) on 
the third shell of a five generation {I — 3, L — 5) Cayley tree with z — 3, and arbitrary 
inner-shell-dependent deposition rate r{l) — OAK The probability reaches a maximum 
of roughly 3.6 x 10"^ at t = 0.037. 

In general, as seen in the previous section for the deposition with constant rates, 
p. ,^ Uie{i,2,..e}Pimi;t) 

P{m,; m^, ..m,; m^; t) = p^g. q. (54) 

where i is the set of interior occupied shells (with at least one empty shell in between), 
and Ltotai is the total number of occupied levels (including the last external shell). 



5. Summary and open questions 

In this paper we presented a class of adsorption models on Cayley trees with constant 
and variable attachment rates, and their possible applications for ionic self-assembly of 
thin films and drug encapsulation. Fig. 9 illustrates a sample comparison between the 
two cases (constant and variable rates) for the probability of having four particles on 
an internal shell £ = 3 for a generation five tree. First note that is set to zero so as 
to block the possibility of dimer deposition (as was considered in Section 3) for better 
comparison with the shell-dependent rate model. The shape of each graph is essentially 
identical. Each has a zero probability at the initial time {at t — the tree is empty), 
slopes up to a maximum, and then completes the trace of a bell skewed slightly to the 
right and approaching zero for large times. Graphs 9(a) and 9(b) differ only in their r2 
value. For Fig. 9(b), r2 = ri = 0.2, meaning we model no electrostatic screening in the 
deposition rates, whereas for Fig. 9(a), r2 = 2ri = 0.4, meaning that a deposition with 
no neighbor is twice as likely as one with neighbors. The third shell of a Cayley tree 
with coordination number z = 3 has twelve open nodes for the m = 4 particles. One 
expects r2 to have a very large impact on the distribution since there are three times 
as many spaces as particles. The maximum probability occurs at t = 0.0565 for (a) 
and t — 1.11 for (b), with a difference of roughly 49.1%. Halving the ri value roughly 
doubles the time scale. The maximum probability remains almost unchanged, with a 
percent difference of about 7.1%. The time scale is what is most affected, with the 
probability of the m = 4 state for the third shell being more spread out for (b). 

The A values in (c) and (d) are chosen to match the r2 values in (a) and (b) 
respectively. This does not make the deposition rate in (c) or (d) equal to those of (a) 
or (b) because r{l) = A', so the rates for the third shell are 0.064 and 0.008 respectively. 
These are much lower and explain the maximum probabilities, which are four and seven 
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orders of magnitude lower than those for the constant rate graphs. The deposition rates 
would be comparable in the first shell, for 1 = 1. For all graphs, we can estimate the 
maximum probability and the time at which it occurs. For variable rates, the peaks 
of the distributions are of the order of 10"'' (for A = 0.4) and, 10~^^ (for A = 0.2), 
practically negligible. The same type of analysis can be done for more general particle 
deposition patterns, where multiple shells are filled with a set number of particles. 

The minimalist cooperative sequential models presented here can be generalized 
to incorporate more realistic elements. In our future work we consider a cooperative 
sequential adsorption model with deposition rates dependent on the number of all 
occupied neighbors. We also want to answer the following basic question: "Under what 
circumstances is a Cayley tree a good enough approximation for a fully two-dimensional 
lattice?" . For the models presented, when ri = rs = ( monomers are blocked to adsorb 
if a neighbor is occupied; no dimers allowed), the particle coverage of the jammed state 
is quite different for the tree (75%) and the lattice (36%). When the occupation of the 
neighboring sites doesn't matter, ri = r2 = 1, the analytical solution for the particle 
density as a function of time matches perfectly the Monte Carlo simulations for a two- 
dimensional lattice. 

For the ISAM problem, the size of the particles has to be addressed, as well 
as possible diffusion of particles on the substrate and evaporation. For the drug 
encapsulation problem, it would be interesting to model the release of the drug molecules 
from the dendrimers and their interactions with the targeted tissue, as well as the 
mobility of the dendrimer branches. These models can also be used to address other 
kinds of problems related to voting behavior and the spread of epidemics. For an 
epidemics model, for example, the "occupied" nodes can play the part of an infected 
individual, while the open nodes are represented by the susceptible individuals. This 
would be the case of a standard susceptible-infected (SI) model. The model can also 
incorporate a recovery rate, to become a susceptible-infected-recovered (SIR) model. A 
version of this model on a Cayley tree is presented in [18j . 
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Figure 1. A Cayley tree of coordination number z = 4 and generation i — Z. 
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Charged Nanoparticles 



Figure 2. Basic steps followed in the layer-by-layer self-assembly process. 
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Figure 3. Sample SEM micrograph of a one-bilayer Si02/ PDDA coating made with 
siHca particles of 45 nm average diameters. 
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Figure 4. (a) Example of a cluster of four connected empty sites, E4; (b) Time 
dependence of the probability of having a cluster of four connected empty sites for a 
Cayley tree with z = 3 and arbitrary deposition rates ri — 0.2, r2 = 0.4 and = 0.1. 
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Figure 5. (a) Time dependence of tlie particle density for a Cayley tree with z — A, 
fi — ^2 — 1, and r3 = 0.1; (b) Steady state coverage for discrete values of 
coordination number z. 
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Figure 6. (a) Time dependence of the probability of having m = 4 particles on the 
inner shell ^ = 3 of a generation-five Cayley tree, with all the other shells empty, for 
z — 4, ri — 0, r2 = 1 and = 0.1; (b) Time dependence of the probability of having 
m = 4 particles on the outer shell L = 5 of a generation-five Cayley tree, with all the 
other shells empty, for z = 4, ri — 0, r2 — 1, and ra = 0.1. 
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Figure 7. Sketch of the drug encapsulation and drug attachment on the outer shell of 
a generation-three dendrimer; different shaped enclosed particles signify the possibility 
of the dendrimer to carry different types of drug molecules. 
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Figure 8. (a) Time dependence of the probability of having m = 2 particles on the 
inner shell ^ = 3 of a generation-five Cayley tree, with all the other shells empty, for 
z = 3 and A = 0.4; (b) Time dependence of the probability of having m = 2 particles 
on the outer shell i = 5 of a generation-five Cayley tree, with all the other shells 
empty, for z = 3 and A = 0.4. 
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Figure 9. Sample comparison of time-dependent probabilities of having m — 4 
particles on the internal shell £ — 3 for a generation-five tree with z = 3 for constant 
rates (a), (b) and variable rates (c), (d). 



